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Human Metapneumovirus Infection Induces Significant 
Changes in Small Noncoding RNA Expression in Airway 
Epithelial Cells 

Junfang Deng^ ^ Ryan N Ptashkin^^ Qingrong Wang\ Guangliang Liu^ ^ Guanping Zhang^ Inhan Lee^ Yong Sun Lee^ and 
Xiaoyong Bao^-^ 

Small noncoding RNAs (sncRNAs), such as microRNAs (mlRNA), virus-derived sncRNAs, and more recently identified tRNA- 
derived RNA fragments, are critical to posttranscriptional control of genes. Upon viral infection, host cells alter their sncRNA 
expression as a defense mechanism, while viruses can circumvent host defenses and promote their own propagation by 
affecting host cellular sncRNA expression or by expressing viral sncRNAs. Therefore, characterizing sncRNA profiles in 
response to viral infection is an important tool for understanding host-virus interaction, and for antiviral strategy development. 
Human metapneumovirus (hl\/IPV), a recently identified pathogen, is a major cause of lower respiratory tract infections in infants 
and children. To investigate whether sncRNAs play a role in hl\/IPV infection, we analyzed the changes in sncRNA profiles of 
airway epithelial cells in response to hl\/IPV infection using ultrahigh-throughput sequencing. Of the cloned sncRNAs, mlRNA 
was dominant in A549 cells, with the percentage of mlRNA increasing in a time-dependent manner after the infection. In addition, 
several hMPV-derived sncRNAs and corresponding ribonucleases for their biogenesis were identified. hl\/IPV M2-2 protein was 
revealed to be a key viral protein regulating mlRNA expression. In summary, this study revealed several novel aspects of hMPV- 
mediated sncRNA expression, providing a new perspective on hl\/IPV-host interactions. 
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Introduction 

Hunnan nnetapneunnovirus (hMPV) is an innportant and 
recently identified nnennber of Parannyxoviridae family that 
also includes respiratory syncytial virus (RSV) and parainflu- 
enza virusJ The genonnic organization of hMPV is analogous 
to RSV, however hMPV lacks the nonstructural genes, NS1 
and NS2, and the hMPV antisense RNA genonne contains 
eight open reading frannes in an order of 3'-N-P-M-F-M2-SH- 
G-L-5'. In ternns of hMPV M2, it encodes two overlapping 
proteins: M2-1 and M2-2. Many clinical studies have shown 
that hMPV is a leading cause of lower respiratory infection 
in infants and children worldwide.^""^ Although nnany stud- 
ies have been carried out, the host-hMPV interaction is still 
poorly understood. 

The discovery that cells produce snnall noncoding RNAs 
(sncRNAs) that regulate gene expression in a sequence- 
specific fashion is a recent breakthrough in biology. 
MicroRNA (nniRNA) is a class of highly conserved sncRNAs 
that usually suppresses gene expression by binding to the 
untranslated regions (UTR) of target nnRNAs, representing 
a new nnechanisnn for the posttranscriptional regulation of 
gene expression.^ The role of sncRNA in regulating antiviral 
innnnune responses has recently ennerged.^'^ Other innportant 



classes of cellular sncRNAs include short interfering RNAs 
(siRNAs), snnall nucleolar RNAs (snoRNAs), and tRNA- 
derived RNA fragnnents (tRFs).^-^° Several viruses, nnajority 
of which are DNA viruses, also produce virus-encoded snnall 
ncRNAs to pronnote their replication or latency.^^'^^ Whether 
the profiles of sncRNAs, both host- and virus-encoded, are 
altered in response to hMPV infection, are unknown. In addi- 
tion, it is not clear whether sncRNAs are functionally innpor- 
tant in controlling cellular responses to hMPV infection. 

High-throughput sequencing technologies, also known as 
deep sequencing or next-generation sequencing technolo- 
gies, deliver fast, accurate, and large volunne sequence data. 
The technology is frequently used for snnall RNA sequencing 
and discovery.^^ We used this advanced technology and bio- 
logical experinnents, such as northern blot and real-tinne poly- 
nnerase chain reaction (PGR), to investigate whether hMPV 
infection results in significant changes in the expression of 
sncRNAs. We found that hMPV infection led to significant 
changes in the expression profile of nniRNAs in hunnan air- 
way epithelial cells. In airway epithelial cells, the changes in 
sncRNA profiles by hMPV infection were different fronn what 
we have recently reported for RSV infection. ^° In response 
to RSV infection, the expression of tRFs is nnost significantly 
affected while there are nnininnal changes in overall nniRNA 
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expression. In contrast to RSV infection, hMPV infection 
did not induce tRFs but gradually enhanced the expression 
percentage of miRNAs in a time-dependent fashion. Expres- 
sion levels of 174 miRNAs were significantly altered over 
a period of 15 hours. Among the altered miRNAs, 142 and 
32 miRNAs were upregulated and downregulated, respec- 
tively, while 72 miRNAs, including miR-30a and miR-1 6, were 
not changed by hMPV infection. Interestingly, even though 
miR-30a and miR-1 6 were not affected by wild-type (WT) 
recombinant hMPV (WT-rhMPV), their expression was sig- 
nificantly induced by a mutant hMPV lacking M2-2 protein 
(AM2-2-rhMPV), demonstrating that the hMPV M2-2 protein 
antagonized the induction of miR-30a and miR-1 6. In addi- 
tion, we identified several hMPV-derived sncRNAs. 

Overall, these data suggest a new way in which paramyxo- 
viruses regulate the host cell response to infection. Given the 
fundamental roles of sncRNAs in regulating the expression of 
genes, including those which are important in cellular signaling 
pathways in other viral infections,^ these findings will provide 
a foundation for future studies investigating the role of various 
miRNAs and its associated pathways in the cellular response 
to this important respiratory virus. The results also potentially 
provide sncRNA-centered information for designing antiviral 
drugs and live attenuated vaccines against hMPV infection. 

Results 

Global sncRNA profile in hMPV-infected A549 cells 

lllumina ultrahigh-throughput sequencing was carried out in 
small RNAs from mock- or hMPV-infected A549 cells to char- 
acterize the global sncRNA expression profile. We obtained 



a total of 58,517,786 sequence reads. These data were pro- 
cessed similarly as described, ^° and is also shown in Figure 1 . 
After stripping adaptor sequences, individual sequences with 
cloning frequency (read numbers) >10 were mapped to the 
human genome (http://hgdownload.cse.ucsc.edu/downloads. 
html) to identify host-derived RNAs and also mapped to hMPV 
genome/antigenome to identify hMPV-derived sncRNAs. 
Sequences mapped to both human and hMPV genomes 
were considered as hMPV sequences when they align better 
to hMPV than to human genome. The small RNAs mapped 
to human genome were further identified as specific human 
sncRNA sequences by comparing them to the miRNA data- 
base (miRBase release 20; http://www.mirbase.org), the rRNA 
database (RDP; http://rdp.cme.msu.edu/), the tRNA database 
(GtRNAdb; http://gtrnadb.ucsc.edu/), the snoRNAs data- 
base (prepared after downloaded sno/miRNA track from the 
UCSC Table Browser), and the Exon-lntron Database (http:// 
www.utoledo.edu/med/depts/bioinfo/database.html). The clon- 
ing frequency of sncRNAs from high-throughput sequencing 
provides a digital measure of their relative expression level. 
To investigate whether hMPV altered the expression profile of 
sncRNAs, we calculated the relative sequencing frequency of 
each sncRNA by dividing its raw read numbers by the total 
read numbers of human genome. 

As shown in Figure 2a, sncRNA composition between 
mock- and hMPV-infected samples is different. The expres- 
sion of miRNAs increased from 66.51% to more than 86% 
after 15 hours of hMPV infection, while small nucleolar RNAs 
(snoRNAs) were gradually altered from 17.91 to 1.49%. 
Minor difference was observed for other types of sncRNAs. 
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Figure 1 Pipeline of analyses of lllumina high-throughput sequencing data. Flowchart of the sequencing data analyses is depicted. 
hMPV, human metapneumovirus; mlRNA, microRNA; rRNA, ribosomal RNA; snoRNA, small nucleolar RNA; tRNA, transfer RNA. 
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Since the percentage of miRNAs clones increased in 
response to hMPV infection, we compare the miRNA profil- 
ing between the mock-infected cells and hMPV-infected cells 
at different time points postinfection As shown in Figure 2b, 
the expression level of 201 miRNAs was upregulated at least 
1.5-fold, while 72 miRNAs was downregulated by <0. 7-fold, 
over a period of 24 hours in infected cells, indicating that 
hMPV induces profound changes in global miRNA expres- 
sion. At 6 hours postinfection, a number of miRNAs were 
already modified in response to the virus. In Supplementary 
Table S1, we have listed miRNAs with early change in their 
expressions and cloning frequency >500. At all time points 
postinfection, upregulated miRNAs were nearly four- to 
fivefold more than the downregulated miRNAs. There were 
also 82 miRNAs that were not affected by hMPV infection 
(Figure 2b). In addition to the change in miRNA expression, 
we also observed 22 hMPV-derived sncRNAs with cloning fre- 
quency >500 in hMPV-infected samples (Table 1).The Inte- 
grative Genomics Viewer was also used to demonstrate the 
distribution of hMPV-derived sncRNAs along the hMPV anti- 
genome (Supplementary Figure S1). These hMPV-derived 
sncRNAs are unlikely to be random cleavage of viral RNAs 
because of the following reasons: (i) these hMPV-derived 
sncRNAs were restrictively antigenome derived, rather 
than genome derived, (ii) not all viral genes were cleaved 



to produce sncRNAs, and (iii) for those viral genes which 
produced sncRNAs, the corresponding regions account for 
a small portion of genes. For example, hMPV polymerase L 
is a gene with more 6,000 nt, however, only few regions have 
their respective RNA fragments. In addition, no quality issue 
was found, as evaluated by 5S rRNA northern hybridization 
(please see the northern blot controls for expression valida- 
tion through the manuscript). 

sncRNA expression validation 

Among the altered miRNA, quantitative real-time PGR 
(qRT-PCR) was used to verify the expression of miRNAs 
which were reported to be affected by other viruses, but with 
functions incompletely characterized or uncharacterized. We 
used U6, a small nuclear RNA, as an internal control for qRT- 
PCR, as its expression was comparable in response to mock 
and hMPV infection (data not shown), although the overall 
expression of snoRNAs was changed by hMPV infection 
(Figure 2a). Let-7f has been reported to regulate viral repli- 
cation or chemokine expression, and its expression is sensi- 
tive to the infection of RSV (a close family member of hMPV) 
and human immunodeficiency virus.^"^"^^ Compared to other 
miRNAs in A549 cells, mock or hMPV-infected, let-7f was the 
second most abundant miRNAs at 15 hours postinfection and 
detectable by northern blot (Supplementary Figure S2). 




Figure 2 Kinetics of changes in sncRNA expression, (a) Human genome-derived small RNAs were sorted according to their origins, and 
their relative abundance was calculated by dividing their frequency numbers by total read numbers, and depicted in pie charts, (b) miRNAs, 
their abundance were changed with a fold change of >1 .5 (upregulated) or <0.7 (downregulated), compared to baseline (uninfected), are 
represented at various times postinfection. Fold change is calculated on relative abundance for each time point compared to the baseline in 
mock infected cells. miRNAs whose expression was not changed are also represented. hMPV, human metapneumovirus; miRNA, microRNA; 
rRNA, ribosomal RNA; snoRNA, small nucleolar RNA; tRNA, transfer RNA. 
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Number of reads Size sequence (nt) IVIost abundant sequence 



Tabie 1 hMPV-encoded sncRNAs with cloning frequency >500 



hIVIPV region hIVIPV gene 

151-183 N 718 

998-1020 N 532 

1977-2007 P 1,330 

2018-2048 P 523 

2261-2294 M 566 

2351-2389 M 747 

2465-2501 M 1,072 

2861-2895 M 1,716 

3158-3189 F 718 

3740-3778 F 2,130 

3801-3829 F 1,132 

4043-4084 F 3,801 

4477-4503 F 728 

4885-4909 M2-1 712 

5090-5129 M2-1 and M2-2 overlapping region 1,537 

5941-5968 SH 741 

6954-6980 Intergenic region between G and L 1 ,214 

7110-7135 L 520 

7135-7163 L 1,248 

7396-7420 L 1 ,026 

7898-7924 L 1 ,226 

13132-13157 L 570 



hMPV, human metapneumovirus; sncRNA, small noncoding RNAs. 

Consistent with our sequencing and northern blot data, 
qRT-PCR assay confirmed that hiVIPV infection significantly 
enhanced let-7f expression (Figure 3a). nniR-452 is inducible 
by influenza H1N1J^ We found its expression was also sig- 
nificantly increased by hMPV (Figure 3a). Whether nniR-452 
plays a role in regulating host response to viral infection is cur- 
rently unknown. We also picked up nniR-374a* and nniR-192 
as the representative of downregulated miRNAs by hMPV for 
expression confirmation. Compared to other downregulated 
miRNAs, miR-192 and miR-374a* were relatively abundant 
with cloning frequency at mock-infected cells above 30,000. 
Both our sequencing data and qRT-PCR revealed that they 
were downregulated by hMPV infection. Similar downregula- 
tion of miR-192 occurs in human hepatoma cells in response 
to hepatitis C virus and hepatitis B virus infection^^'^^ and 
the involvement of miR-192 in enhancing replication of the 
hepatitis C virus has also been shown. Of note is the fact 
that the information on the regulation of miR-374a* by viral 
infection is very limited except their induction by hepatitis C 
virus.^^ Our data demonstrate that miR-374a* was reduced 
by hMPV infection (Figure 3a). In this study, we also used 
qRT-PCR to confirm the expression of miR-30a and miR16, 
whose expression was not affected by hMPV infection. 

As shown in Table 1, hMPV infection resulted in the 
expression of hMPV-derived sncRNAs as well. It is com- 
monly believed that negative-sense RNA paramyxoviruses 
do not produce viral miRNAs, as their replication does not 
occur in the nuclear compartment, where miRNA biogen- 
esis is initialized. ^^'^^ However, several other sncRNAs are 
known to be generated by alternative pathways, and some 
RNA viruses produce virus-derived sncRNAs.^'^'^'^^ To confirm 
the sequencing data of hMPV-derived sncRNAs, we used 
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northern blot to detect the two hMPV-derived sncRNAs. As 
shown in Figure 3b, two sncRNAs derived from the P and 
L, two genes of hMPV RNA-dependent RNA polymerase 
complex, were detectable by northern blot. The band corre- 
sponding to 201 8-2048 was slightly higher than 71 35-71 63, 
which is consistent with their sequencing data showing 
that the overall length of 2018-2048 was longer than that 
of 7135-7163. In addition, we used 6- andl 0-day film expo- 
sure to acquire the bands corresponding to the regions of 
7135-7163 and 2018-2048, respectively, which again is in a 
line with the sequencing data suggesting that there was more 
induced 7135-7163 than 2018-2048. 

Function of let-7f in hIVIPV replication 

Among the differentially expressed miRNAs, we further 
explored let-7f function. Similar to hMPV infection, RSV 
infection has been shown to enhance the expression of let- 
7f expression as well, and the enhancement actually favors 
RSV replication.^"^ To determine whether hMPV-induced let-7f 
is functional, we first investigated whether hMPV replication 
is sensitive to the treatment of let-7f inhibitor. Commercially 
available miRNA inhibitors are usually chemically modified, 
single-stranded nucleic acids designed to specifically bind to 
and inhibit endogenous miRNA molecules. These ready-to- 
use inhibitors can be introduced into cells using transfection 
or electroporation parameters similar to those used for siR- 
NAs and enable detailed study of miRNA biological effects. 
A549 cells were transfected with 150 nmol/l let-7f inhibitor 
or its control for 24 hours, followed by mock or hMPV infec- 
tion at multiplicity of infection of 2. After 1 hour, the cells 
were replaced with fresh media and incubated for additional 
15 hours. Cells were then scraped, followed by sonication 



Molecular Therapy — Nucleic Acids 



hMPV Regulates sncRNA Expression 

Deng etal 



1.6 r 



1.4 



^ 1.0 



0.8 



0.6 



Let-7f 

miR452 

miRSOa 

miR16 

miR374a* 

miR192 




2018-2048 



hMPV (hours): 0 6 15 0 6 15 0 6 15 0 6 15 0 6 15 0 6 15 



5s rRNA 




15 hours 



Figure 3 Experimental validation of sncRNA expression, (a) miRNA validation. Total RNA from indicated treatments in A549 cells was 
harvested at indicated time points postinfection, followed by reverse transcription using aTaqMan microRNA reverse transcription kit (Life 
Technologies). Quantitative real-time polymerase chain reaction was used to validate the expression of representative miRNAs. U6 snRNA was 
used as an internal control. All reactions were done in a 1 0-|jl reaction volume in triplicate. Data are representative of two to three independent 
experiments. *P< 0.05, **p< o.OI , relative to the first bar in each group, (b). hMPV-derived sncRNA expression. A549 cells, mock- or hMPV- 
infected, were harvested at 15 hours postinfection, followed by total RNA extraction and loading to a denaturing polyacrylamide gel for northern 
hybridization using indicated probes. 5S rRNA is shown for equal loading. Data are representative of two to three independent experiments. 
hMPV, human metapneumovirus; miRNA, microRNA; rRNA, ribosomal RNA; snoRNA, small nucleolar RNA; tRNA, transfer RNA. 

and centrifugation to release the intracellular viral particles. 
The infectious hMPV particles were nneasured by innnnunos- 
taining as we previously described and expressed as plaque- 
fornning unit/nni (the nunnber of infective particles within 1 nnl 
of the sannple).^^-^^ We found that cells transfected with let-7f 
inhibitor had a fold change of 6 in hMPV replication connpared 
to control oligotransfected cells (Figure 4a). A real-tinne PGR 
for hMPV N gene copy nunnbers also showed sinnilar trends 
(data not shown). nniRNA nninnics are snnall, chennically nnodi- 
fied double-stranded RNAs that nninnic endogenous nniRNAs 
and enable nniRNA functional analysis by upregulation of 
nniRNA activity. In contrast to the effect of let-7f inhibitor on 
hMPV replication, lef-7f nninnic significantly decreased hMPV 
replication 3.3-fold (Figure 4b).Transfection of let-7f-specific 
nninnic or inhibitor did not result in nnore cytotoxicity in cells 
than control oligos as the trypan blue staining did not reveal 
the difference in dead cell nunnbers between let-7f-specific 
and control oligo-treated cells (data not shown). These data 
suggest that hMPV-nnodulated let-7f has an antiviral role. 

The role of hMPV M2-2 protein in miRNA expression 

M2-2 is an hMPV protein that regulates viral RNA 
synthesis.2^-2° Recently, we have discovered that hMPV M2-2 
is a potent virulent factor that blocks nnitochondrial antiviral 
signalling-nnediated antiviral signaling leading to inhibited 
activity of two transcription factors: NF-kB and interferon reg- 
ulatory factor (IRF)-3, and, subsequently reducing production 
of type I interferon (IFN).^^ Since nniRNAs are RNA molecules 
and its regulation has been shown to be regulated by tran- 
scription factors, such as NF-kB and IRF-3, and IFNs,^^"^"^ 
we investigated whether M2-2 regulates the synthesis of 
nniRNAs. In brief, A549 cells were infected with reconnbi- 
nant hMPV (rhMPV), either WT or a nnutant lacking M2-2 
(AM2-2). At 15 hours postinfection, cells were harvested for 
RNA extraction followed by nniRNA quantification by qRT- 
PCR. Interestingly, we did not observe any effect of M2-2 
on the expression of let-7f and nniR-452, two hMPV-induced 
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Figure 4 Let-7f inhibits human metapneumovirus (hl\/IPV) 
replication, (a) Knockdown of let-7f inhibits hMPV replication. Let- 
7f inhibitor was transfected into A549 cells using lipofectamine 2000. 
At 24 hours post-transfection, cells were mock infected or infected 
with hMPV at multiplicity of infection (MOI) of 2. Cells were harvested 
at 6 or15 hours postinfection, and viral infectious particles were 
measured by immunostaining. (b) Let-7f mimic were transfected 
into A549 cell. At 2 hours post-transfection, cells were mock infected 
or infected with hMPV at MOI of 2. At 6 or 15 hours postinfection, 
viral titration was assayed. Data are representative of two to three 
independent experiments. *P< 0.05, relative to the white bar at each 
time point postinfection as indicated, pfu, plaque-forming unit. 
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Figure 5 M2-2 regulates microRNA (mlRNA) expression, (a) A549 cells were transfected with plasmid expressing M2-2 or its control 
(0.1 |jg/well). At 30 hours post-transfection, the cells were mock infected or infected with rhMPV, wild type (WT), or AM2-2, at multiplicity of 
infection (MO!) of 2. At 15 hours postinfection, total RNA was prepared followed by quantitative real-time polymerase chain reaction (qRT-PCR) 
to measure the expression of miR-30a and miR-16. (b) U4A cells were infected with rhMPV at MO! of 2 for 15 hours, followed by total RNA 
preparation and qRT-PCR to measure the expression of miR-30a and miR-1 6. U6 snRNA was used as an equal loading internal control. Data 
are representative of two to three independent experiments. *P< 0.05, **p< o.OI , relative to the first bar in each group, (c) The total protein 
from U4A cells was also subjected to Western blot to assess the viral gene expression. hMPV, human metapneumovirus. 



nnlRNAs. M2-2 also did not affect the expression of nniR-192 
and nniR-374a*, two hMPV-decreased nnlRNAs (data not 
shown). As nnentioned above, M2-2 deletion significantly 
enhances IFN-a/p secretion and the activation of NF-kB and 
IRF-3 in response to hMPV infection.^^ Connparable expres- 
sion of let-7f, nniR-452, nniR-192, and nniR-374a* in WT- or 
AM2-2-infected cells suggests that the effect of hMPV on the 
expression of these nniRNAs does not depend on IFN signal- 
ing or the activation of NF-kB and IRF-3. 

Notably, connpared to WT infection, AM2-2 infection 
resulted in an increased expression of nniR-16 and nniR-30a 
in airway epithelial cells, although the expression of both 
nniRNAs was not affected by the infection of naive hMPV or 
WT-rhMPV (Figures 3a and 5a). To confirnn the role of M2-2 
in nnodulation of nniR-16 and nniR-30a, the M2-2 protein was 
exogenously expressed. We found that M2-2 expression alle- 
viated enhanced nniR-16 and nniR-30a expression following 
rhMPV-AM2-2 infection (Figure 5a), dennonstrating a role of 
M2-2 in regulation of these two nniRNAs. We also found that 
the inhibition of M2-2 on nniR-16 induction was innpaired in 
U4A cells, a cell line lacking IFN signaling because of JAK-1 
knockout, ^^'^^ suggesting that M2-2-nnediated nniR-16 expres- 
sion depends on IFN signaling (Figure 5b). Connparable 
nniR-16 expression in WT- and AM2-2-infected U4A cells also 
suggested that IRF-3 and NF-kB are not innportant for nniR-1 6 
induction, as M2-2 deletion resulted in nnore activated IRF-3 
and NF-kB in U4A cells, akin to the regulatory role of M2-2 
in the activation of IRF-3 and NF-kB in airway epithelial cells 
(data not shown). In contrast to nniR-16, nniR-30a seenned to 
be IFN independent. In the absence of JAK-1 , AM2-2 induced 
nnore nniR-30a than WT virus in U4A cells (Figure 5b). 
Besides direct inhibition of antiviral signaling pathway by 
M2-2, M2-2 also suppresses other viral gene expression to 
facilitate evasion of hMPV fronn host innnnunity.^^ Therefore, it 
is possible that M2-2 regulates nniR-30a through nnodulating 
other viral gene expression. However, our results show that 
M2-2 deletion does not change hMPV protein expression in 
U4A cells (Figure 5c), suggesting a direct regulatory function 
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Figure 6 Human metapneumovirus (hl\/IPV)-derived sncRNA is 
generated by Dicer and RXN1. (a) A549 cells were transfected 
with 120 nmol/l of sIRNA against indicated proteins or control 
siRNA as a negative control. At 48 hours post-transfection, the cells 
were mock- or hMPV-infected for 24 hours. Total RNAs were then 
subjected to northern hybridization as described in Figure 3b; 5S 
rRNA expression is shown for equal loading, (b) The suppression of 
target mRNAs by each siRNA was confirmed by quantitative real- 
time polymerase chain reaction measurement. The values are the 
relative expression levels of each mRNA, with the value of control 
siRNA set as one. Data are the representative of two independent 
experiments. 

of M2-2 in nniR-30a expression. We show, for the first tinne, 
that M2-2 regulates nnlRNA expression. 

Biogenesis of hIVIPV-derived sncRNAs 

Most studies on the roles of virus-derived sncRNA in control- 
ling viral replication were perfornned with double-stranded 
DNA viruses. It is connnnonly believed that negative-sense 
parannyxovi ruses do not produce viral nniRNAs as their 
replication does not occur in the nuclear connpartnnent. 



IVIolecular Therapy — Nucleic Acids 



hMPV Regulates sncRNA Expression 

Deng etal 



where miRNA biogenesis is initialized. However, several 
sncRNAs are known to be generated by alternative pathways 
and some of them exert miRNA-like functions.^ Indeed, our 
deep sequencing data demonstrate that there were sev- 
eral hMPV-encoded sncRNAs from hMPV infections. These 
sncRNAs could be viral miRNAs generated via noncanoni- 
cal pathways, or other forms of sncRNAs, such as siRNA or 
recently identified novel microRNA (miRNA)-like molecules 
generated from virus and host tRNA transcripts. To define 
their functions in the future, it is important to identify the 
enzymes responsible for their biogenesis. As mentioned, 
hMPV is a single-negative-stranded RNA virus. Therefore, 
the transcription and replication of hMPV likely happens 
in the cytoplasm, and hMPV-derived sncRNAs cannot be 
generated by nuclear RNases, such as Drosha and RNase 
P, two enzymes critical for miRNA biogenesis and tRNA 
cleavage, respectively. ^^"^^ As expected, we did not observe 
the changes in the expression of hMPV-derived sncRNAs 
by Drosha silencing (Figure 6), reinforcing our view that 
the biogenesis of hMPV-derived sncRNAs was distinct 
from that of siRNA/miRNA. To identify host protein(s) 
responsible for hMPV-derived sncRNA generation, we paid 
close attention to the nucleases whose expression was 
significantly induced by hMPV, using previously collected 
microarray data.^^ One particular enzyme called exoribonu- 
clease 1 (XRN1) was identified as a candidate. Since Dicer 
is a cytosolic nuclease in the RNase III family that cleaves 
double-stranded RNA and pre-miRNA into short double- 
stranded RNA fragments, and we observed that its gene 
transcription is enhanced by hMPV infection, it was also 
identified as a candidate. We found that both XRNI and 
Dicer are necessary for the expression of hMPV-derived 
sncRNA upon hMPV infection (Figure 6a). The efficiency 
of knockdown was confirmed by qRT-PCR measurement 
of each mRNA (Figure 6b) and western blot (data not 
shown). Collectively, our data suggest that the production 
of hMPV-derived sncRNAs upon hMPV infection specifi- 
cally required XRNI and Dicer. 

Discussion 

The discovery of posttranscriptional control of genes by 
sncRNAs is a key advance in biology. sncRNAs, such as 
miRNA and recently identified tRFs, are gene-silencing 
mediators, which have an important role in the regulation of 
many important cellular process, including cellular differenti- 
ation, innate immunity, apoptosis, and oncogenic transforma- 
^jQP, 9,10,42 During host-virus interactions, many viruses have 
been reported to change miRNA expression profile to modu- 
late host responses to viral infection. ^^-^^ Whether sncRNAs 
are changed by hMPV infection is currently unknown. In addi- 
tion, there is limited knowledge on the relationship between 
paramyxovirus-induced host responses and sncRNAs, 
either host or virus encoded. In this study, we found that the 
expression of miRNAs among sncRNAs in A549 cells was 
increased post-hMPV infection. In addition, the changes in 
cellular sncRNA profiling in response to hMPV are different 
from that of RSV Following RSV infection, tRFs are the major 
type of sncRNAs significantly affected, while the change in 



miRNAs was minimal. ^° However, in the context of hMPV 
infection, hMPV infection did not induce tRF expression but 
changed miRNA profiling significantly. A total of 375 known 
miRNAs were cloned in A549 cells, among which, 142 and 
32 miRNAs were upregulated and downregulated by hMPV 
at 15 hours postinfection, respectively. 

As mentioned, there are very limited studies on the inter- 
action between miRNAs and paramyxoviruses. Hendra 
virus, a highly pathogenic zoonotic paramyxovirus, has been 
reported to make use of miR-146a to negatively regulate 
NF-kB activity via targeting ring finger protein 11, to favor 
its replication.^^ RSV downregulates miR-221 , which leads to 
enhanced expression of the nerve growth factor in human air- 
way cells, and subsequently favors viral replication by interfer- 
ing with the apoptotic death of infected cells.^^ RSV also uses 
let-7f induction to promote its replication.^"^ Since let-7f was 
the second most abundant miRNAs with cloning frequency 
>1 00,000 reads in uninfected cells and was also significantly 
induced by hMPV infection, we investigated whether let-7f is 
functionally important. In contrast to the stimulus role of let-7f 
in RSV replication, we found that hMPV-induced let-7f inhib- 
ited hMPV replication (Figure 4). Similar inhibitory effect of 
let-7f on viral replication has been shown for several other 
viruses, ^^"^^ suggesting that let-7f is an important and com- 
mon regulator of viral replication. In addition, in these experi- 
ments, the lactate dehydrogenase assay did not reveal the 
cell toxicity by transfection reagent (data not shown). The 
opposite regulatory effect of let-7f on hMPV and RSV repli- 
cation suggests that the role of miRNA in regulating viral rep- 
lication is pathogen dependent. There are several possible 
reasons why miRNA functions differently in response to dif- 
ferent viral infections. First, miRNA may target different host 
genes in the context of hMPV and RSV infection. Second, 
since the sncRNA-viral RNA interaction has been reported 
to be important in viral RNA synthesis and stability, ^° the dif- 
ferent role of let-7f in RSV and hMPV replication may result 
from the difference in the binding of let-7f to RSV and hMPV 
RNA. To investigate whether miRNAs are able to interact with 
hMPV genome/antigenome, we mapped miRNAs to the viral 
genome/antigenome.^^ Interestingly, we found that the first 
13 bases of let-7f match with hMPV genome 9871-9883, 
where RNA-dependent RNA polymerase L is located. How- 
ever, such interaction did not exist between let-7f and RSV 
genome.^^ The mechanism underlying the effect of let-7f on 
hMPV will be studied in the future. In the future, we will also 
investigate the functions of other miRNAs whose expression 
was affected by hMPV infection. 

In this study, we also demonstrated that hMPV M2-2 pro- 
tein inhibited the expression of miR-30a and miR-16, as 
AM2-2 infection resulted in a significant induction of miR-30a 
and miR-16, while WT infection did not result in induction 
(Figure 5a). By target prediction algorithms (TargetScan and 
miRanda), several targets were hit. Given the fact that miR-1 6 
and miR-30a were not sensitive to naive/WT hMPV infection, 
we selected those targets, whose expression is not changed 
by naive hMPV in our previous microarray gene transcrip- 
tion data, as targeting candidates.^^ For example, we chose 
cysteinyl-tRNA synthetase (CARS), KIAA1949, CDC37L1, 
and runt-related transcription factor 1 (RUNX1) as target 
candidates of miR-30a, because their expression remained 
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unaffected by hMPV infection. RUNX1 actually is required 
for human immunodeficiency virus infection. Whether these 
targets favor hMPV infection and whether hMPV uses M2-2 
to suppress these miRNAs to prevent downregulation of their 
targets will be studied in the future. It is a common trend 
that miRNAs with significant change in the expression by 
viral infection are always paid close attention. This example 
reminds us that we may need to study those miRNAs that 
are unchanged by viral infection, as they are probably sup- 
pressed by certain viral components to favor viral infection, 
which, in this case, is hMPV M2-2 protein. 

We also identified several hMPV-encoded sncRNAs and 
corresponding ribonucleases for their biogenesis (Figure 6). 
Whether these sncRNAs play a role in mediating viral replica- 
tion and innate immune responses need to be experimentally 
validated. If so, such virus-encoded small RNAs would be ideal 
tools for controlling viral replication, as they are nonimmuno- 
genic. They may also regulate multiple host transcripts to vary- 
ing degrees. Therefore, silencing the action of virus-encoded 
sncRNAs may enable host cells or the immune system to gain 
control over the virus and even to eliminate it. Due to the scant 
homology between viral sncRNA and host cell miRNAs,^^ anti- 
sncRNA oligonucleotides might have minimal off-target effects 
and great therapeutic potential. A recombinant virus system for 
hMPV is well established in our laboratory Any modification to 
the recombinant genome, such as modifying sequences cor- 
responding to virus-encoded sncRNAs, should help develop 
live attenuated vaccines against hMPV infections. In addition, 
hMPV, as a negative sense RNA virus, might be a suitable 
model to study the biogenesis of virus-encoded sncRNAs for 
those viruses that do not replicate in the nuclear compartment. 

In summary, deep sequencing results provided us an over- 
all big picture without any priori assumptions and guided us 
which sncRNAs species to be explored further, although data 
bias may exist. Within reasonable variance, we confirmed the 
deep sequencing results with qRT-PCR and northern blot 
experiments which overcame the data bias problems. Over- 
all, our ultrahigh-throughput sequencing study provides novel 
information on the regulation of sncRNAs by hMPV infection 
and can serve as a foundation for future investigations on the 
role of the identified miRNAs and hMPV-derived sncRNAs 
in initiating antiviral as well as innate and adaptive immune 
responses to this important pathogen. In addition, sncRNAs 
identified to be important for hMPV replication have poten- 
tials to serve as therapeutic targets, suggesting a potential 
translational application for this study. 

Materials and methods 

Cell lines and virus. A549, human alveolar type ll-like epi- 
thelial cells, and MK-2 cell line (both from ATCC, Manas- 
sas, VA) were maintained as described. ^^'^^ U4A cells, a 
JAK-1 -deficient fibrosarcoma cell line (kindly provided by 
Dr George Stark, Cleveland Clinic Foundation, Cleveland, 
OH), were cultured as previously described.^^'^^ hMPV, 
either naive or recombinant derived from the Canadian iso- 
late hMPV83 (CAN 97-83), was grown in MK-2 cells, puri- 
fied by sucrose-gradient and titrated by immunostaining in 

MK-2 cells using homemade rabbit anti-hMPV antibody as 
described. 26-28,54-56 



Viral infection. A549 cells were infected with hMPV, at mul- 
tiplicity of infection of 3, in serum-free medium containing 
antibiotics and 1 |jg trypsin/ml for all experiments, unless oth- 
erwise stated. Cells treated with infection media plus sucrose 
were defined as mock infected. Mock- or hMPV-infected 
cells were scrapped into the medium, followed by centrifu- 
gation and RNA extraction using TRIzol Reagents according 
to manufacturer's instruction (Invitrogen, Carlsbad, CA) as 
described. ^° 

Deep sequencing and RNA mapping. A549 cells after mock 
or hMPV (CN 97-83) infection were harvested for total RNA 
preparation. The RNAs were delivered to Eureka Genomics 
(Houston, TX) for small RNAs isolation, directional adap- 
tor ligation, cDNA library construction, and sequencing 
using a Genome Analyzer llx (lllumina, San Diego, CA). 
About 2,103 Mb of sequence data with a total of 58,517,786 
sequence reads were generated for mock- and hMPV-infected 
samples, using 36b single-end sequencing reads. The cDNA 
sequences were aligned to both human (hg19) and hMPV 
(CAN97-8) reference genomes using NovoalignCS. The 
alignments were carried out with options to strip 3' adapter 
sequences, for a minimum of 15 good quality bases per read, 
and a tolerance of one mismatch. After initial alignment, fur- 
ther processing was performed using in-house programs and 
SAMtools.^° First, uniquely aligned reads and sequences 
aligned to more than one genome location (ambiguously 
aligned reads) were separated. The ambiguously aligned 
reads were then randomly assigned to one location and com- 
bined with uniquely aligned reads. Sequences with reads >1 0 
were used for the downstream analysis shown in Figure 1. 
We assigned the sequences mapped to both hgl 9 and hMPV 
to only one of them with better BLAST score not to duplicate 
the read number. We performed several normalization meth- 
ods for comparison including total reads of human genome, 
total number of each species (such as total miRNA reads or 
total tRNA reads), and upper quartile normalization. While 
detailed fold changes differ, the overall trends are similar. 
Considering the shifting of sncRNAs composition with hMPV 
infection, we used values normalized by total human reads 
in the results. 

RNA interference. Transfections of siRNA into A549 cells 
were carried out at a final concentration of 120 nmol/l, tar- 
geting Drosha, Dicer, 5'-3' exoribonuclease 1 (XNR1), or a 
scrambled negative control (Sigma, Houston, TX), using lipo- 
fectamine 2000 (Invitrogen) according to the manufacturer's 
recommendations. Fourty-eight hours later, A549 cells were 
mock infected or infected with hMPV for 18 hours at a multi- 
plicity of infection of 3. 

Expression confirmation of sncRNA. Northern hybridiza- 
tion for sncRNAs were performed as described. ^'^^ Briefly, 
RNA was separated in 15% denaturing polyacrylamide gel 
with 7mol/l urea and then transferred to a positively charged 
nylon membrane (Amersham Biosciences, Piscataway, NJ). 
The membrane was hybridized with 32P-labled probes in 
ULTRAhyb-Oligo solution (Life Technologies, Grand Island, 
NY), followed by washing according to the manufacturer's 
instruction. 
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qRT-PCR was also used to validate mlRNA expression data 
from high-throughput sequencing. Extracted total RNA was 
reverse transcribed into single-stranded cDNA with a Taq- 
Man microRNA reverse transcription kit (Life Technologies). 
Real-time PGR was performed using first-strand cDNA with 
TaqMan Fast Universal PGR Master Mix. The assay numbers 
for the miRNA endogenous control (U6 snRNA) and miR- 
NAs were as follows: U6 snRNA (001 973), miR-1 6 (000391 ), 
let-7f (00382), miR-1 92 (000491), miR-1 94 (000493), miR- 
374a*(002125), miR-30a (000417), and miR-452 (003229). 
Quantitative PGR was performed on an Applied Biosystems 
Prism 7500 Fast Sequence Detection System (Applied Bio- 
systems). Quantitative PGR parameters for cycling were as 
follows: 95 °G for 20 seconds, 40 cycles of PGR at 95 °G for 
30 seconds and 60 °G for 30 seconds. All reactions were 
done in a 10-|jl reaction volume in triplicate. The mRNA and 
miRNA expression levels were determined using the 2~A^'^ 
method. 

Statistical analysis. Statistical significance was analyzed 
using analysis of variance. P value of less than 0.05 was con- 
sidered significant. Mean ± standard error is shown. 

Supplementary material 

Figure S1. Distribution of hMPV-derived sncRNAs. 
Figure S2. Let-7f validation. 

Table S1. miRNA with total sequencing reads >500, and fold 
induction >1.5 or fold reduction <0.7 by hMPV at 6 hours 
postinfection. 
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